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AN ANALYSIS OF THE FAC SIMILE -CAMERA RESPONSE 


TO RADIANT POINT SOURCES 

By Friedrich O. Huck, Stephen J. Katzberg, Daniel J, Jobson, 
and Carl L. Fales, Jr. 

Langley Research Center 

SUMMARY 

In addition to imaging the surrounding terrain, planetary lander cameras may also 
be used to survey the stars to aid in locating the lander site. The response of the facsim- 
ile camera, which was selected for the Viking lander missions to Mars, to a radiant point 
source is formulated and shown to result in a statistical rather than deterministic signal. 
The signal statistics are derived and magnitudes are evaluated for the brighter visual and 
red stars. The probability of detecting the resultant statistical signals in photosensor 
and preamplifier noise and the associated probability of false alarms are also determined. 

INTRODUCTION 

Although the primary purpose of the Viking lander imaging experiments on Mars is 
to visually characterize the landing sites, it will also be valuable to observe the brighter 
stars and planets. These observations could improve the accuracy of locating the landers, 
and, thus, aid in correlating Viking orbiter and lander science data (ref. 1). 

A logical first step in preparation for such an experiment is to determine the 
response of the facsimile camera, which was selected for the Viking landers, to radiant 
point sources. At a first glance it may appear that this problem simply reduces to the 
detection of deterministic signals masked by noise. Actually, however, as will be shown, 
the facsimile camera generates a signal which is statistical rather than deterministic. 

This paper formulates and evaluates the facsimile-camera response to radiant point 
sources, the statistics of the resultant signals, and the probability of detecting the brighter 
visual and red stars in photosensor noise. 

SYMBOLS 

Ac camera objective lens aperture area, cm2 

AAc elemental area of lens, cm2 



H hypothesis 

spectral irradiance from point source, W/cm^-jnm 

normalized spectral irradiance from blackbody radiator or Planckian emitter, 
pm-1 

I signal current, A 

i integer 

spectral radiant intensity, W/sr-pm 
K magnitude of photosensor signal current, A 

K' normalized magnitude of photosensor signal current, A 

k peak spectral irradiance, pW/cm2 

L lens -aperture point spread function 

N number of picture elements (pixels) per scan line 

n noise 

P probability 

P^ spectral radiant power, W/pm 

AP^ elemental spectral radiant power, W/pm 

p probability density 

spectral responsivity of photosensor, A/W 
S photosensor-aperture point spread function 

s signal 
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T integration period of photosensor current, sec 

Tg sampling period of photosensor current, sec 

t time, sec 

X angular vertical sampling interval, rad 

X normalized period during which signal current is integrated (see eq. 6(c)) 

Y angular azimuth samplir^ interval, rad 

0 ! = 1 - |3 (see eq. 6(a)) 

jS photosensor integration period normalized to the sampling interval X (see 

eq. 6(b)) 

A magnitude of sampled signal, A-sec 

6 delta or impulse function 

? probability of detecting a star 

7] normalized signal current threshold, A 

6{x) Heaviside step function, 1 for x > 0, and 0 for x = 0 

X wavelength, pm 

^ probability of false alarm 

p angular radius of photosensor aperture, rad 

a root-mean- square (rms) value of electronic noise current, A 

T optical transmittance 

X vertical angle of lens coordinates, rad (see fig. 2) 


3 



Xg vertical angle of lens-mirror scanning coordinates, rad (see fig. 2) 

Ax vertical width of photosensor-aperture section, rad (see fig. 8) 

angular azimuth stepping interval, rad 
azimuth angle of lens coordinate, rad (see fig. 2) 
x{/g azimuth angle of lens-mirror scanning coordinates, rad (see fig. 2) 

Ai// azimuth width of photosensor-aperture section (see fig. 8) 

w angular distance between photosensor aperture and lens optical axis, rad 

n(x) rectangular function, 1 for | X | = ^, and 0 elsewhere 

oo 

ni(x) sampling or comb function, ^ 6(x - i) 

i=-oo 

Subscripts: 

a atmosphere 

c camera 

n noise 

o photosensor aperture 

0 no star sample 

1 one star sample 

2 two star samples 

1,2 at least one star sample 

A bar over a symbol represents its average value. 

A dot over a symbol represents the rate of motion. 
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ASSUMPTIONS 


The analysis presented in this paper is subject to the following assumptions: 

(1) The lens point spread function is small compared to the photosensor aperture 
and can be neglected. The constraints imposed by this assumption on the results of this 
paper, especially for the Viking lander camera, are discussed in appendix A. 

(2) The photosensor -aperture shape is circular. Other photosensor -aperture 
shapes, which may often be preferable (ref. 2), would generate different signal statistics. 

(3) A silicon photosensor is used. 

(4) Only one point source, or star, with sufficient irradiance to be detectable is 
present within a selectable field of view, or set of image data. 

IMAGING PROCESS 

A basic configuration of the facsimile camera is illustrated in figure 1. Radiation 
from the object field is reflected by the scannii^ mirror, captured by the objective lens, 
and projected onto a plane which contains a photosensor aperture. The photosensor con- 
verts the radiation falling on the aperture into an electrical signal which is then amplified 
and sampled for digital transmission. As the mirror rotates, the imaged object field 
moves past the aperture, thus permitting the aperture to scan vertical strips. The cam- 
era rotates in small steps between each vertical line scan until the entire object field of 
interest is scanned. (In addition, the Viking facsimile camera has the capability of 
repeating vertical line scans in a fixed azimuth position.) 

Two coordinate systems must be accounted for as illustrated in figure 2. One coor- 
dinate system with an origin at the center of the mirror (fig. 2(a)) accounts for the cam- 
era scanning geometry, and the other coordinate system with an origin at the center of the 
lens (fig. 2(b)) accounts for the imaging geometry. The vertical angle through which the 
mirror has scanned the object field is labeled Xg ” where Xg is the mirror scan- 
ning rate, and the azimuth angle through which the camera has rotated is labeled . 

O 

The axis of azimuth rotation is also the optical axis of the objective lens and is labeled z 
in figure 2. The location of an object point (with the mirror optically unfolded) and the 
corresponding image point is measured as an angle from this common axis at the center 
of the objective lens. One angle, measured alor^ the direction in which the imaged object 
field is moved by the scanning mirror, is labeled x> the other angle, measured nor- 
mal to this direction, is labeled i//. The angular location of the photosensor aperture is 
given by the angles Xq ’^o' should perhaps be emphasized for clarity that the 

angles Xg (tig. 2(a)) and x (fig- 2(b)) occur in the same plane, but that the angles 
(fig. 2(a)) and xp (fig. 2(b)) occur in planes which are normal to each other. Since the 


5 




Figxire 1.- Basic facsimile -earner a configuration. 

photosensor aperture is always located at or near the optical axis, only small angles 
are of interest; thus, (tan x^tan \J/) = (x,^)- 

Let J;y be the radiant intensity of a point soufee, and let 5(x,^) be the delta 
(ref. 3) or impulse (ref. 4) function. Then, as shown in figure 3 , the elemental spectral 
radiant power AP;^^ incident upon an elemental area of the camera lens AAq is 


= “2 "^aX^^c cos 
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Figiore 3.- Defining diagram for emitted and captured radiation. 
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where is the transmittance of the optical path between source and lens and 

oj w ^ ;^2. Since the image of a point source must be near the optical axis in order 
to pass over the photosensor aperture (i.e., w ~ 0), and since the distance L of interest 
is great compared to the lens-aperture dimensions, AA^. may be replaced by the lens- 
aperture area A^.. Also, by substituting the spectral irradiance of the point source 

/where for its radiant intensity J^, the foregoing equation becomes 


- ^cHx^aX 


where AP^ has been replaced by P^, the total spectral radiant power incident on the 
camera lens. 

The transfer of the spectral radiant power of the point source which is captured by 
the objective lens P^ S(x,V/) to the photosensor-aperture plane may be expressed by the 
convolution 



6(x' dx' 


= Px 


where L{x,4^) is the lens point spread function. The effect of the scanning mirror is to 
shift the imaged radiant power over the photosensor-aperture plane, and the effect of the 
camera rotation between line scans is to advance the imaged radiant power in integral 
steps normal to the shifting direction. The process by which the photosensor converts 
this radiant power into an electrical signal may, therefore, be expressed by 


/ \ ^00 



where is the transmissivity of the optical path in the camera, is the respon- 

sivity of the photosensor, and S(x,i/^) is the angular extent of the photosensor aperture. 

The symbol rn(;///Y) is the sampling (ref. 3) or comb (ref. 4) function. This func- 
tion is essentially an infinite sum of delta functions whose spacings in this case are 
related to the angular azimuth stepping interval by Y = ^tg cos Xg if Xg is mea- 
sured from a plane normal to the optical axis of the objective lens. 


8 



The photosensor -aperture shape is assumed to be circular and can be expressed as 


Six,^) = 1 



Six,^) = 0 


(Elsewhere) 


where p is the angular radius of the aperture and is small. 

The lens point spread function, or diffraction pattern, is assumed to be small com- 
pared to the photosensor aperture and can, therefore, be approximated by the impulse 
function; that is, L(x,J/^) ~ Equation (1) becomes then 


l(Xo,K,Xs^,^) = Jq Px^cX^X ^ ni(|] = ^ ™ 


( 2 ) 


where K is the magnitude of the signal current given by 


pOO pOO 

K = H^V^CX^X <S) 

Equation (2) may be visualized with the help of figure 4 as presenting an infinite 
number of "tracks” of a point image, with spacing Y between successive tracks. These 
tracks generate a signal current ^ duration that a track crosses 

the photosensor aperture; otherwise, ^(Xo’^o’^s^’'^) ~ 


X 



Figure ij;.- Illustration of point-source image tracks in the photosensor-aperture plane. 
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The signal current is passed to a running mean integrator (ref. 5) for sampling - 
a process which may be formulated as follows: Let T be the period of signal integra- 
tion so that XgT becomes the effective angular interval over which integration occurs; 
and, similarly, let Tg(>T) be the period of sampling so that X = XgTg becomes the 
effective vertical sampling interval. The sampled current may then be expressed as 





-i- ml -A 

XY lx Y/ 


(4a) 


or, by using the rectangular function 



it can be expressed as 



(Elsewhere) 







(4b) 


where * represents the convolution operation. 

Similarly, as before, this expression may be visualized with the help of figure 5 as 
presenting again an infinite number of tracks of a point image with spacing Y. However, 
although a signal current with magnitude K is generated for the duration that a track 
crosses the photosensor aperture, this current is integrated only during the interval XgT, 
once per sampling interval X. Consequently, the value of the sampled signal depends not 
only on the magnitude of the current generated by a point source but also on where and 
when the star image crosses the photosensor aperture, thus giving rise to a signal which 
is statistical rather than deterministic in nature. 
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X 



Figure 5.- Illustration of point-source image tracks and effective signal- integration 
intervals in the photosensor-aperture plane. 


SIGNAL MAGNITUDE AND STATISTICS 

The foregoing analysis has shown that the signal generated by a facsimile camera, 
when scanning a radiant point source, depends not only on the source irradiance and cam- 
era sensitivity but also on the camera scanning geometry, photosensor -aperture shape, 
and signal processing characteristics. This section determines the magnitude of the 
signal current generated by the brighter stars and the statistics of the sampled signal 
current. 


Signal Magnitude 

Reference 6 has shown that stars behave essentially as Planckian emitters (black- 
body radiators), and that, therefore, the spectral irradiance for a star may be obtained 
from its visible radiance and effective temperature. In order to generalize results, it is 
convenient to rewrite equation (3) in the form 


K' = = k r r[Rx dX (5) 

TcAc 

where K' is the signal current normalized with regard to optical design parameters of 
the facsimile camera (optical transmittance Tq and lens -aperture area Ag), is 
the normalized blackbody radiation, and k the peak spectral irradiance for a star. It is 
assumed that the atmosphere transmittance is unity, as is approximately true for 

the Martian atmosphere within the spectral responsivity of the silicon photosensor (except 
near the terminator) , and that the camera transmittance is not wavelength dependent. 
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The spectral responsivity of a typical silicon photosensor is plotted in figure 6, and the 
resultant photosensor current K' is plotted in figure 7. A list of the brightest visual 
and red stars for which the peak spectral irradiance k was calculated in reference 6 is 
given in table I, together with visual magnitude, effective temperature, and normalized 
photodiode current K'. 

TABLE L- CHARACTERISTICS OF THE BRIGHTEST VISUAL AND RED STARS AND 
NORMALIZED SILICON PHOTOSENSOR SIGNAL CURRENT 


Stars 

Visual 

magnitude 

(a) 

Effective 
temperature, K 

(a) 

Peak spectral 
irradiance, k, 
pW/cm2 

(a) 

Photosensor 
current, K', pA 

Sirius 

-1.60 

11 

200 

40 

1.21 

Canopus 

-.82 

6 

200 

2.9 

.28 

Vega 

.14 

11 

200 

8 

.24 

Capella 

.21 

4 

700 

2.4 

.30 

Arcturus 

.24 

3 

750 

2.8 

.35 

Rigel 

.34 

13 

000 

9 

.18 

Achernar 

.60 

15 

000 

20 

.27 

Beta Centauri 

.86 

23 

000 

30 

.10 

Altair 

.89 

7 

500 

2.2 

.16 

Betelgeuse (variable) 

.92 

2 

810 

3.8 

.36 

Pollux 

1.21 

3 

750 

1.2 

.15 

Antares 

1.22 

2 

900 

2.3 

.23 

Alpha Crucis 

1.61 

2 

810 

3.2 

.30 

Mira (variable) 

1.70 

2 

390 

3.8 

.26 

Beta Gruis 

2.24 

2 

810 

1 

.10 

R. Hydrae (variable) 

3.60 

2 

250 

.8 

1 

.05 


^Data taken from reference 6. 


Signal Statistics 

As is illustrated in figure 8, often only part of the signal current K is actually 
integrated and contributes to the sampled signal magnitude. Figure 8(a) depicts the photo- 
sensor aperture divided into rectangular sections of width A\p and length and fig- 

ure 8(b) shows two of these sections in more detail. The path of the star image across 
each one of the rectangular aperture sections is depicted by a thin line, and an arbitrarily 
selected integration period y„T is depicted by a heavy line superimposed on each thin 
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line. It can be seen that either two, one, or no signals may be generated and that a sig- 
nal may vary in magnitude from essentially zero to K^Xgl^X^ The number of star sig- 
nals obtained and their normalized magnitude (K = 1) depend on two independent events; 
namely, (a) where (along the azimuth direction) the star image passes across the photo- 
sensor aperture with respect to the aperture center, and (b) when (along the vertical scan 
direction) the star image passes across the photosensor aperture with respect to the time 
of signal current integration. 

Within a single line scan, the signal current statistics of interest are the probability 
of generating two, one, or no signal, and the cumulative probability distribution, probabil- 
ity density, and average value of the integrated magnitude of the sum of two signals and of 
at least one signal. These statistics are presented here without derivation (analytical 
details are given in appendix B) and are based on the following assumptions and 
constraints: 

(1) A star image may cross the photosensor aperture anywhere along azimuth with 
equal probability. Only one crossing is considered here; multiple crossings are consid- 
ered in the next section. 

(2) A star image may cross the photosensor aperture anytime with respect to the 
start of integration with equal probability. 

(3) The angular vertical sampling interval X is assumed to be equal to the 
photosensor-aperture diameter 2p, as is generally the case. 

In order to simplify the notation, it is convenient to make the following four 
definitions: 

(1) The normalized period during which integration does not occur is given as 


X - XgT XgT 

a = ^ = 1 - 

X X 


(6a) 


(2) The normalized period during which integration does occur is given as 




X 


= 1 - 


a 


(6b) 


(3) The normalized period during which integration occurs while the image of a point 
source is passing across the photosensor aperture is given as 


X = 



C0 


where 0 < c = 1; hence, 0 < x = ^. 
14 


(6c) 




(4) A sampled signal normalized with respect to X is defined by the symbol A; 
if two signals are to be distinguished, they are noted as Aj and A 2 . The magnitude of 
a signal may vary, therefore, from zero to jSK; that is, 0 < A/K = jS. 

Probability of signal generation .- The probabilities of generating two, one, or no 
star signals are, respectively (see appendix B), 





■^Q! - Q!\l - a^') 


(7a) 


Pi(a>o} = |/| + 2|3 + 2q!\ 1 - 0-2 - 2 cos" Iff) 


Pq{a=o} = -[20! + cos"^Q! - - - a\ 
2 \ 2 


1 - 


(7b) 


(7c) 


The probability of generating at least one star signal can be obtained from the previously 
given results as 


Pi,2(a>0} = P2{Ai>0,A2>0} + Pi(a>0} = 1 - Po(a=o} (7d) 


Figure 9 illustrates the variation of these probabilities with the normalized integration 
period jS. 



Figure 9.- Variation of the probability that two, one, at least one, or no signal 
sample is generated with normalized integration period p. 
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Cumulative probability distribution . - Two probability distributions are of interest 
(see appendix B) and are given as follows; 

(1) The probability that the sum of two star samples, Aj^ + A 2 , is larger than xK 
is given as 


|x<^A i+A 2 ^ j K=/^ = i cos"l(a + x) - (a - x)^l - (a + x)^ 


(8a) 


P2|^^^<(Ai+A2)/Kj = 0 


(8b) 


(2) The probability that at least one star signal, Aj or A 2 or both, is larger 
than xK is given as 


Pj 2(x<VKSp> =1 


cos'^x - cos" t(a + 2x) + (2j3 - 3x)^l - x^ + (a + 2x)^l - (a + 2x)^ 


Pi^2{x<A/KS^}=|. 


cos- lx + (2/8 - 3x) yi - x2 


(x < 8/2) 


(8/2 S X £ 8) 


> (8c) 


Pi^2(/3Sx<a} = 0 


(8d) 


When X < /3/2, it is possible for signal Aj or A 2 or both to be larger than xK; when 
|3/2 = X = |3 it is possible only for either signal Ai or A 2 (but not both) to be larger 
than xK. 

Figure 10 represents these probability distributions for several normalized integra- 
tion periods /3 and unity signal current K = 1. It may be noted that equation (8a) prop- 
erly reduces to equation (7a) for x = 0, thus giving the probability of generating two star 
signals; and, similarly, equation (8c) reduces to equation (7d) for x = 0, thus givit^ the 
probability of generating at least one star signal. 

Probability density .- The probability-density function of a continuous random vari- 
able, in this case x, is defined as the derivative of the cumulative probability distribu- 
tion. (Because the distributions given by eqs. (8) increase with decreasing x, the nega- 
tive value of the derivative must be used here; i.e., p(x) = -dP(x)/dx.) Consequently, the 
density functions correspondii^ to equations (8) are given as follows; 
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(a) Probability that the sum of two star samples is larger than x. 



(b) Probability that at least one star sample is larger than x. 


Figure 10.- Cumulative probability distribution of star samples for several values 

of normalized integration periods 3- 

(1) The sum of two star samples: 

flO . X) (9a) 

\jl - (o + x)^ 

e(x) = 0 if X S 0. 


P2|x;^Ai+A2)y 


K) = 


where 0(x) = 1 if x > 0 and 
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(2) At least one star sample: 




1 - {a + 2xY 


(x<^/2) 


1- 


X 


> (9b) 


1 - x" 


P^^2 (x;A/K) = — ■ 00 - X) + -|(cos-l^ - ^\ll - 6(x - (3) 0/2 S x S 1) 


J 


Figure 11 presents these probability densities for several normalized integration 
periods jS and unity signal current K = 1. It may be noted that the impulse function 
represents the probability that a star signal current is integrated over the complete inte- 
gration interval 0 § x = 1. 

Average value .- The average, or expected, value of a random occurrence, in this 
case the probability density function p(x), is defined as 

J OO 

xp(x) dx 

_oo 


Consequently, the average values of the occurrences described by equations (9) are given 
as follows: 

(1) The sum of two star samples: 



+ A2^ = ot sin'^Q! + -i^2 + 1 - q!^ - — a 


(10a) 


(2) At least one star sample: 

(10b) 

Figure 12 presents the variation of these average values with the normalized inte- 
gration period /3; the signal current is assumed to be unity, that is, K = 1. 
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0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 


X 

(b) At least one star sample. 

Figure 11.- Probability-density function of star samples for several values 
of normalized integration period 3. 
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Figiare 12.- Variation of average star-sample magnitude plotted against 

integration period for K = 1. 

DETECTION AND FALSE ALARM 

Besides converting radiant power into an electrical signal current and amplifying it, 
the camera photosensor and preamplifier also generate a noise current. This section 
determines the probability of detecting star signals which are characterized by the sta- 
tistics formulated in the foregoing section and are masked by this noise. 

General 

Based on the foregoing results, the imagery data contained in a single line scan may 
be searched for two signals, one signal, or at least one signal of a single star. The latter 
approach has the highest probability of success and is, therefore, pursued here. For sim- 
plicity, analytical results are evaluated only for signal integration intervals equal to one- 
half of the sampling interval (i.e., for a = |3 = 0.5), as is the case for the Viking lander 
cameras. 

By using a standard statistical approach (see, for example, ref. 7), the following 
simple hypotheses may be set up for each signal sample or picture element (pixel): 

Horn 

Hj : s + n 

In words, hypothesis Hq stipulates that only noise is present, and hypothesis Hj stip- 
ulates that a star signal is also present. 
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The noise generated by the photosensor and preamplifier is random and can be 
expressed by the Gaussian density function (ref. 5) as 


Pn(x) 


1 -x^/2a2 

- ' — e ' 

\2-na 


( 11 ) 


where a is the root-mean-square (rms) value of the noise current (in amperes). 

The probability density for signal plus noise p., « (x) is given by the convolution 
of the signal and noise probability density functions and can be written as 




( 12 ) 


where pj^ 2 (^) 1® given by equation 9(b) and is plotted in figure 11(b). The two probabil- 
ity density functions Pj^(x) and Pj^ ^ ^^(x) for hypotheses Hq and Hj, respectively, 
are illustrated in figure 13. 


CO 


S3 

0 > 


■8 


CXi 


.4 


.3 


.2 


.1 



X 


Figure 13 •- Probability density of noise and signal plus noise 
for a signal-to-noise ratio of k/o = 6. Also Illustrated 
is the probability of star detection ^ and false alarm | 
for a threshold t]. Normalized signal- integration period 
of p = 0.5. 
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Let it be assumed now that for any signal current x 
threshold value r}{i.e., x = t]), as illustrated in figure 13, 
X > T], Hj is accepted. The probability of accepting 
given by 


smaller than or equal to some 
Hq is accepted, and that for 
when Hi is, in fact, true is 



Similarly, the probability of accepting when Hq 


is true is given by 


(13) 



The probability often called the power of the test, is referred to in this paper as 
the probability of detecting a star; the probability often called the level or size of the 
test, is referred to as the probability of false alarm. 

Since signal plus noise are digitized for transmission as n binary digit (n-bit) 
words, the threshold value rj is forced to be one of 2^^ discrete levels. In the Viking 
lander cameras, for example, n = 6 and the spacing between quantum levels is selected 
by gain settings. 

It should be noted that 


pOO pOO 

J Pl,2,n(x) dx = ] Pi ^ 


where is the probability that at least one star signal will be generated, and 

Pq {a= 0^ is the probability that no star signal will be generated; both probabilities are 
given by equations (7). Consequently, the probability that a star which is present will not 
be detected is P^ i “ C ~ 

No single approach for searching for a star signal in a set of image data may prove 
consistently superior to another. Nevertheless, it appears to be generally advantageous 
to obtain several successive line scans across a star, display this set of line scans con- 
tiguously and carefully registered, and then search for a line pattern across the line scans. 
Successive line scans across a star may be obtained by allowing the facsimile-camera 
azimuth rotation to overtake slowly the apparent star motion (due to planetary rotation) , 
or vice versa. The Viking lander cameras also have a so-called repeated line-scan mode 
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which inhibits azimuth rotation without otherwise charging the camera imaging process, 
thus providing another opportunity to obtain several line scans across a star. If suc- 
cessive line scans across a star can be obtained, then the approximate slope of the line 
pattern may possibly be estimated from considerations of apparent star motion with 
respect to the camera scanning rate and geometry. If this approach is pursued, then the 
most useful a priori information may well be the probability of detecting at least one star 
signal and at least one false alarm in each one of the line scans. 

Since it was assumed that only one star with sufficient irradiance to be detectable 
is present with a selectable field of view, or set of image data, the probability of detecting 
at least one star signal in a single line scan is given by equation (13). This probability is 
plotted in figure 14 for several threshold ratios 77 /a against the signal-to-noise ratio 
K/a. 



Figure ill-.- Probability of detecting at least one star signal -within 
a single line scan plotted against slgnal-to-nolse ratio for se-v- 
eral thresholds r\/a. Normalized signal-integration period of 

p = 0.5. 
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A false alarm, however, can occur at each pixel in the line scan. These occur- 
rences are independent events and may be recognized as an example of Bernoulli trials 
(ref. 8). Thus, the probability of detecting i false alarms in N pixels is given by 




and of detecting at least one false alarm is given by 


N 


^ = 1 - (1 - 

i=l 


(15) 


( 16 ) 


where ^ is given by equation (14). The latter probability is plotted in figure 15 for the 
same threshold ratios r]/o as used in figure 14. 



Figure 15 •- Probability of at least one false alarm in W pixels 
for several thresholds Tj/a. 
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Application 


In order to estimate the probability of detecting one of the brighter stars, it is 
necessary only to look up the normalized signal current K' from table I, calculate the 
actual signal current K = as given by equation (5), measure the camera photo- 

sensor and preamplifier rms noise current a, and enter the signal-to-noise ratio K/a 
into figure 14. Consider the Viking facsimile-camera design (ref. 1). The objective 
lens-aperture diameter is 0.95 cm and the total optical transmissivity is 0.7, so that 
AcTc = 0.5 cm^. Hence, from table I, the signal current, for example, for the star 
Canopus becomes K = (0.5)(0.28 pA) = 0.14 pA. The total noise current of the photo- 
sensor and preamplifier may be taken as 4 fAy/(|Hz (ref. 9), so that a = 0.02 pA for 
the 25 Hz bandwidth of the camera slow-scan mode. Hence, the signal-to-noise ratio is 
K/a = 0.14/0.02 = 7. Entering this ratio into figure 14 yields a probability of detection of 
about 0.7 for a threshold of 77/a = 2. 

In order to estimate the probability of false alarm from figure 15, it is necessary 
only to determine the number of pixels per line scan N (i.e., vertical field of view 
divided by instantaneous field of view) which may be expected to contain a star signal 
from astronomical considerations. For example, for N = 8 pixels per line scan, the 
probability of at least one false alarm is almost 0.2 for a threshold of 77/0 = 2. The 
effect of data transmission errors on false alarm is formulated in appendix C. 

CONCLUDING REMARKS 

The facsimile -camera response to a radiant point source, such as a star, was for- 
mulated and shown to result in a statistical rather than deterministic signal. The signal 
statistics were derived by using several simplifying assumptions, of which the most con- 
straining one is that the lens point spread function is small compared to the photosensor 
aperture and can be neglected. Probable signal magnitudes were determined for the 
brighter visual and red stars, by assuming that the facsimile camera uses a silicon photo- 
sensor. Since these signals are masked by photosensor and preamplifier noise, their 
probability of detection and the associated probability of false alarms were also formu- 
lated and evaluated. Curves are presented which depict these probabilities for several 
threshold signal-to-noise ratios. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 15, 1973. 
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APPENDIX A 


LIGHT DISTRIBUTION OF A POINT-SOURCE 
IMAGE NEAR FOCUS 


This paper assumes that the image of a point source is essentially a delta or 
impulse function in the photosensor -aperture plane. This assumption cannot strictly be 
true, of course, because of diffraction at the objective lens, and it is further violated 
when defocus blur is present. The purpose of this appendix is threefold: first, to illus- 
trate the actual size of an error-free diffraction image of a point source near the lens 
focal point; second, to suggest an approach which would account for an extended image of 
a point source; and third, to discuss the effect of the assumption of a point-source image 
on the results presented in this paper. 

The light distribution near focus in an error -free diffraction image is given in ref- 
erence 10. For the present purpose it may be sufficient to consider the total illumina- 
tion of a point-source image which would fall within a circular aperture of radius v 
centered on the lens optical axis at a distance u from the lens focal point. Figure 16 
presents a contour line for a somewhat arbitrarily selected radius v which would con- 
tain 80 percent of the total illumination. The normalized defocus distance is given by 
(see refs. 10 and 11) 


u 



(Al) 


where D is the lens -aperture diameter, I is the photosensor -aperture distance from 
the lens, and AZ is the photosensor -aperture distance from geometric focus. For com- 
parison with the radius v, the normalized photosensor -aperture radius of the facsimile 
camera is given by 


V 


c “ 



(A2) 


where R^ is the actual photosensor-aperture radius. 

Consider now the design parameters of the Viking lander camera (ref. 1) listed as 
follows: 
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■ APPENDIX A - Continued 

Lens focal length, f, cm 5.25 

Lens -aperture diameter, D, cm 0.95 

Photosensor -aperture radius, Rg, 19.0 

Photosensor-aperture distance from lens, I, cm^ 5.27 

Photosensor-aperture distance from geometric focus for a point source located 
at infinity, AZ, cmt> 0.02 


^For the high -resolution photosensors; the low-resolution photosensor apertures 
are three times wider. 

^For the high-resolution photosensor closest to the lens focal point. 

Substituting the foregoing parameters into equations (Al) and (A2) yields v^. = 15 and 
u = 15 for the average wavelength of the silicon photosensor, X = 0.7 jum. 



Figiire l6.- Contour line for 80 percent of total illumination of a point-source 
image which falls within circles of radius v centered on the lens optical 
axis at distance u from the lens focal point (from ref. 10) . 
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APPENDIX A - Concluded 


Figure 16 shows that v = 3 for u = 0 and that v is, therefore, less than one- 
fifth of the photosensor -aperture radius Vc- Thus, the assumption of a point-source 
image should yield fairly accurate results for this imaging situation. 

As can also be seen from figure 16, v = 14 for u = 15 so that v is then 
approximately equal to Thus, the assumption of a point-source image is no longer 

proper, and some of the results presented in this paper are not valid for this imaging 
situation. 

One approach to account for defocus blur would be to convolve the defocused lens 
point spread function (given, for example, in refs. 10 and 11) with the photosensor- 
aperture point spread function as formulated by equation (1), and to use the resulting 
function as the effective (and enlarged) photosensor -aperture point spread function. The 
image of a point source may then again be regarded as a delta function, and the signal 
statistics may be derived as outlined in the paper and detailed in appendix B. 

Upon some reflection it may be visualized that the effectively "enlarged" photosen- 
sor aperture would now allow more than two signal samples of a single point source to be 
generated within a single line scan. The signal statistics presented in figures 9 to 12 
would, therefore, significantly change — but with one important exception. This excep- 
tion pertains to the probabilities associated with the generation and detection of at least 
one star signal. For the Viking lander camera, for which the normalized integration 
period /S is equal to 0.5 (as used for figs. 13 and 14), the probability of generating at 
least one star signal (see fig. 9) is 0.98. This probability could not, of course, be signifi- 
cantly improved upon. Consequently, the final result presented in figure 14 — the prob- 
ability of detecting at least one star signal - is not expected to change significantly in the 
presence of some defocus blur. The results presented in figure 15 - the probability of 
at least one false alarm - do not depend on any of these considerations. 
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APPENDIX B 


DERIVATION OF SAMPLED-SIGNAL-CURRENT 
STATISTICS FOR POINT OBJECTS 


The purpose of this appendix is to present in detail the derivation of some of the 
statistics which are used in this paper to describe the sampled signal current generated 
by point objects. The number of star signals obtained and their magnitude depend, as 
previously described in the paper and illustrated in figure 8, on two independent events: 

(1) where (along the azimuth direction) the star signal passes across the aperture center, 
and (2) when (along the vertical scan direction) the star image passes across the photo- 
sensor aperture with respect to the time of integration. The signal statistics derived 
here are (a) the probability of generating two, one, or no signal sample, and (b) the cumu- 
lative probability distribution of the integrated and sampled-signal magnitudes. The lat- 
ter statistics, (b), may be more clearly defined as the probability that at least one star 
sample (or the sum of two star samples) is larger than some value x, where x is an 
independent variable. It follows from this definition that (b) reduces to (a) for x = 0. 
However, since the purpose of this appendix is primarily to show how the sampled-signal 
statistics were derived, it is convenient to derive first (a) and then (b). 

All assumptions and constraints are stated in the main part of this paper, together 
with illustrations of the results. Two exceptions are made for this appendix: (1) Instead 
of requiring that the effective vertical angular sampling interval X be equal to the cir- 
cular photosensor diameter 2p, it is allowed that X = 2p; and (2) the signal current K 
is assumed to be unity for convenience of notation. 


Probability of Generating Two, One, or No Signal Sample 

Two situations may be recognized from figure 8: (1) When Ay- > X - XoT, then 

1 s 

either two or one star samples occur: and (2) when Ay- = X - XaT, then either one or no 
star sample occurs. For any rectangular section i for which Ay^ > X - x„T, the prob- 

1 '’S 

ability of obtaining two star signals is given by 


, . 1 f^Xi+XsT-X 1 


- (X - XsT 


and of obtaining one star signal is given by 

^X 




^Xi+X< 


dx=^ 

T-X X 


(2X - XgT) - AXi 


(Bl) 


(B2) 
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APPENDIX B - Continued 


Similarly, for any section i for which Ay. = X - Xc,T, the probability of obtaining one 

A '’S 

star signal is given by 

= S (B3) 

and of obtaining no star signal is given by 

4I^x.Av. ^ 

Now, let 1 S i = n - 1 be the range for which Ay. > X - x„T, and let n = i = m 

be the range for which Ax^ = X - x^T. The summation of all rectangular sections over 

s 

the complete photosensor aperture may then be expressed as 



It follows from the foregoing definition for n that 
Ay„ = X - Xc.T 

An As 

and from the geometry of the circular photosensor aperture that 
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APPENDIX B - Continued 


Equating these two equations for AXj^ and solving for n yields 


n 


1 

2 



X - 


(B8) 


It follows similarly from the foregoing definition for m that 




= 0 


and from the geometry of the circular aperture that 


Ax 


m 




Again, equating the two equations for Ax^^^ and solving for m yields 



P 

Axl/ 


(B9) 


In the limit as A\j/ approaches zero, the summations in the probability equa- 
tions (B5) to (B7) become integrals and are given as follows: 

For generating two star samples, 


P2(Ai>0,A2>0} 


Ai 

Xp 



where 


and 



if/ s lim 4^ . = lim 4/ = 

A r\ 11 J. A . 1 . r\ 


Ai^-0 


Aj ^-0 




32 



APPENDIX B - Continued 


Solving this equation yields for a ' = 1 


P2(Ai>0,A2>0} 



■la' - 



where 


(BIO) 


a 


X-XsT 

2p 


2p 


a 


By this relation, a' may exceed unity if X > 2p since 0 ^ a S 1. Thus far it has 
been implicitly assumed that ;// = 0 or a'S!. If 0!'>1, that is, if X - XgT > 2p, 
two star signals cannot be generated and, hence. 


P2(o!’>i} ^ 0 

Similarly, for generating one star sample for o;’ S 1, 


(Bll) 




where \ ^ given previously. Solving this equation yields for a = 1 


Pi{a>o> =£ 


-2 cos'^a' - 2 q!'[i - '\|l 


ce'2 .X + 5 

/ P 2 


(B12) 


where a' is given previously. Also, for a’ > 1 


‘ t(-“' * i " f) 


(B13) 


Finally, for generating no star sample for o; ^ 1, 




X-x,T)(p-^)-j!2|x|d^ 
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APPENDIX B - Continued 


Pq{a=0} =^^os"1q!' - 


a 


,2 


. 2«' - I 


and for a' > 1, 



- Xs'^ - 2jx|)di^ = 



(B14) 


(B15) 


Cumulative Probability Distribution of Signal Samples 

Three cumulative probability distributions are derived: (a) the sum of two star 
signals when two star signals are present, (b) one star signal when two star signals are 
present, and (c) one star signal when one star signal is present. The distribution for one 
star signal when either one or two star signals are present (i.e., of at least one star sig- 
nal) is the sum of cases (b) and (c). 

Three situations may be recognized similarly as before: 

(1) When AXj. > X - XgT + xX, then either two or one star samples definitely occur 
with a total magnitude A (where A = Aj + A 2 for two samples) greater than x, 

(2) When xX < A\, = X - v„T + xX, only one star sample may occur with a signal 
magnitude greater than x. 

(3) When xX = Ax^, no star sample can occur with a signal magnitude greater than 

X. 

The condition of two star samples in case (1) may be separated into two distinct 
parts: (la) When Ay. > X - x^T + 2xX and two star samples occur, a portion of the 
interval X has both signals Aj and A 2 each greater than x; and (lb) when 
X - x„T + xX < Ax, = X - XcT + 2xX, then either one star sample or the other may occur 

o X S 

with a signal greater than x but not both. 

From figure 8, for any rectangular section i for which Ax, > X - x_T + xX, the 

X '’S 

probability that there are two star samples and that the magnitude of their sum Aj^ + A 2 
is larger than some value x is given by 


P 2 ,{x<Ai+A2S/3'> = 


dx-i 

^ X 


iXi - (x - X , t ) 


(B16) 


and 


P 2 i{(J’fe<Ai+A2> s 0 


(B17) 
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APPENDIX B - Continued 


where 


^' = /3 - 1 + 


X 


(B18) 


and for = 0, two samples cannot exist. The probability that there are two star sam- 
ples present, at least one of which is greater than x, is given by 


^Xi+XgT-X 

= i dx = i(AXi ^ - X) 


where > X - + 2xX, and 

'"X s 


/ Axj-xX-X.x,T 

P2_i{x<ASp'}=-U dxH-\ 


AXi+XgT-X 


dx 


(B19) 


P2^i{x<AS^’}=|(AXi + XsT-X-xx) 

where X - v^T + xX < Ay. = X - y„T + 2xX, and finally 

''D X D 


(B20) 


P2,i{/3’&<A> = 0 


(B21) 


The probability that only one star sample is present and that the signal magnitude is 
larger than x is given by 

P^ i{x<AS^’> = i r dy = 1(2X - y T - Ay.) (B22) 

XJAy.+XgT-X XV V 


where Ay. > X - y T + xX, and 
X s 


P, . {x<AS^'> - if dx = i(Ax, . X,T - 2xx) 


(B23) 


where xX < Ay^ = X - XgT + xX, and 

Pf i(^'Sx<A} = 0 


(B24) 
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APPENDIX B — Continued 


All rectangular sections may again be summed over the photosensor aperture. 
The new limits of summation may also be obtained similarly. Let ^3 

denote the boundaries between the aforementioned regions (la) and (lb), (lb) and (2), 
and (2) and (3), respectively. Thus, 




1^/2 = P\1 - (a’ +x^ 


1//3 = p\|l - X 


.t 2 


V 

where x' =~x. 

2p 

Furthermore, when 1 S a' + 2x' and 1 > a’ + x' ^or equivalently ^ = a + 2x 
and ^ when 1 = 0 !'+ x', then xj/^ = O; and when 1 = x’ 

(or 2p = xX), then ^ ^ 

The probability that there are two star samples and that the magnitude of their sum 
Ai + A 2 is larger than some value x is given by 


P^{x<Ai+A2^^'} 

P2(x<Ai+A2=/3'} 

P2{x<Ai+A2S^'} 



and for all values of a' and x’, 


(a' +x' < 1) 


(a’ +x’ 

<1) 

(B25) 

(a' + x' 

^1) 

(B26) 


P2(|3’gx<Ai+A2} = 0 (B27) 

The probability that there are two star samples present, at least one of which is 
greater than x, is given by 
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APPENDIX B - Continued 




( 2 |x| + X^T -x)d^. 2(21^1 , - X - xx) 


(a' + 2x' < 1) 


P^{x<A^^'} =^{a +2x')^l - (a' + 2x')^ - 2(a' + x’) U ~ + x’)^ 


cos"^(q!' + 2x’) + 2 cos"^(q!' + x') 


(a' + 2x' < 1) (B28) 


P 2 {x<A^^'} = ± 2[2|x| + XgT - X - xX 


d;^ (a' + 2x' SI; 0 ;' + x' < 1) 


P2(x<AS3') 


XL 


cos"^(a!' + x') - {a' +x')\ 


1 - {a' +x')2 


(a' + 2x' S 1; a' + x' < 1) (B29) 


P2 (x<A^| 3'} = 0 


(a’+x’Sl) (B30) 


and, finally, for all values of a' and x'. 


P2(^'Sx<a} h 0 


(B31) 


The probability that only one star sample is present and that the signal magnitude 
is larger than x is given by 


Pl{x<AS3'} 


(2X - %T - 2|x|) dV. H. (2|x| + X3T - 2xx’ 


^3 


di// 


{a' + x' < 1) 
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APPENDIX B - Continued 


P.{x<A^|3’} =1 


X 


2{a' + x')^l - (o!' + x')^ + ^-3x' - 2a' + 


.,2 


2 cos“^(q!' + x') + cos"^x’ 


(q!'+x'<1) (B32) 


Pl{x<ASr> = A r® (2|x| + - 2xx) a*- 


P^(x<AS|3'} = 
Pj{x<Ag/3'} = 0 


_p_ 

cos~^x' + 

(-3x' - 2 q!' 

1 - x’2 

X 





(q!' + x' = 1; x' < 1) 

{a' + X' S 1; X’ < 1) (B33) 

(x' ^ 1) (B34) 


and, finally, for all values of a' and x', 

P^(^'^<a} = 0 (B35) 


At the beginning of this appendix it was noted that X = 2p in order that the deriva- 
tions have some generality. However, in the main part of the paper it was assumed that 
X = 2p. Thus, 


0!' — Q! 


)3' ^ /3 

x' — X 


and the pertinent equations of (BIO) to (B35) become 


P2{Aj> 0,A2>0^ =i^cos"^Q! - 0!|l - 


Pl(A>0}=|| + 2^ + 2a^l-a2 


1 _ qi 2 _ 2 C0s“^Q!l 


(B36a) 

(B36b) 


38 



APPENDIX B - Concluded 


Pq{a=o} = ^(20! + cos~^a! - - Q!^l - q ;2 


P2(x<Ai+A2^^} 


- 1 0 
cos (o! + x) - (a - x)\ 1 - (Q! + x) 


P2{i3^<Ai+A2} = 0 


P2(x<ASp> = 1 


{a + 2x)\l - {a + 2x)^ - 2(o! + x)\|l - (o; + x)^ 


cos"^(q! + 2x) + 2 cos~^(a! + x) 


(x < /3/2) 


P2(x<A=^} 


cos ^(fl! + x) - (q! + x)\l - (a + x) 


(^/2 g X < |3) 


P2(/3gx<A} = 0 


Pj{x<ASp}=i 


2(o! + x)\]l - {a + x)^ + (2^ - 3x)\|l - x2 - 2 cos“^(o! + x) + cos"^x 


-1 


(B36c) 

(B37a) 

(B37b) 

(B38a) 

(B38b) 

(B39a) 


Pj{|3gx<A} = 0 


(B39b) 



APPENDIX C 


TRANSMISSION ERRORS 

In addition to the camera photosensor and preamplifier noise, the imagery data are 
also susceptible to transmission errors. The purpose of this appendix is to formulate 
the effect of these errors on the probability of generating false alarms. 

In general, a transmission error can chaise the value of an n-bit pixel word which 
is below the threshold value t] so that it is still received below this threshold (in which 
case the decision that no star is present is unchanged), or so that it is received above the 
threshold (in which case a false alarm is generated). Similarly, a transmission error 
can change the value of a word which is above the threshold value rj so that it is still 
received above this threshold (in which case the decision that a star is present is 
unchained) , or so that it is received below this threshold (in which case a star signal is 
lost). In addition, it may sometimes be possible to predict that a star signal is very 
unlikely to exceed a determinable upper level. Some of the most significant bits of a 
n-bit word can then be disregarded, thus reducing the probability of false alarm. 

Only the probability of false alarm for low transmission-error rates will be con- 
sidered here. The probability of false alarm is higher than the probability of losing a 
star signal simply because a line-scan data record contains many pixels of which only 
one or two represent a star signal. 

Let 6 be the probability that a binary digit (bit) is in error. Then, the probability 
that at least one bit of a n-bit word (or the n least significant bits of a larger word) is 
transmitted incorrectly is given by 

y = 1 - (1 - e)"^ (Cla) 

\ 

Generally, e < 0.01 and n < 10 so that equation (Cla) can be approximated as 

y = ne (C lb) 

In words, this approximation assumes that it is very unlikely that more than one trans- 
mission error occurs in a single word. 

Now, let it be assumed that a n-bit word (or the n least significant bits of a larger 
word) which has a magnitude below or equal to the threshold value 17 is transmitted. 

The probability that this word will then be received with a magnitude above 77 is given by 
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APPENDIX C - Concluded 


^t = en 

= e^n - 3.32 log^^ rj) (C2) 

The probability that the word transmitted is, in fact, below the threshold value 77 
is 1 - 4, where | is the probability of false alarm given by equation (14). Since this 
probability is independent of the probability the probability that a false alarm will 
occur because of a transmission error is given by 

r = 4t(l - (C3) 

Finally, since a false alarm can result independently either because of camera 
noise (with probability or because of a transmission error (with probability |'), the 
probability of false alarm in any given pixel becomes 

i” = r + 4 

r = e(n - 3.32 log^Q t7)(1 - ^ (C4) 

Also, the probability of at least one false alarm in N pixels becomes 

^i’,N = 1 - (1 - (C5) 
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